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Semldlrect  methods  utilize  direct  (non-iterative)  linear  solvers  to  solve 
nonlinear  problems  in  a non- time-like  iteration.  The  semldlrect  methods  are 
applied  to  various  fluid  dynamics  problems.  In  each  case,  the  direct  linear 
solver  used  is  some  variation  of  a marching  method.  ^ 
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STATEMENT  OF  THE  PROBLEM  STUDIED 

The  scope  of  the  contract  work  was  to  develop  and  Improve  efficient 
semi direct  methods  for  solving  the  discretized  nonlinear  partial  differen- 
tial equations  of  various  fluid  dynamics  problems.  The  basic  concept  of 
"semidirect''  methods  Is  to  use  direct  (non-iterative)  fast  solvers  for 
linear  equations  to  solve  nonlinear  problems  using  non-time-like  Iterations. 
The  direct  solver  used  herein  Is  some  variation  of  marching  methods;  the 
research  progress  thus  involved  both  the  nonlinear  semidirect  methods  and 
the  tool  which  these  require,  the  direct  marching  methods  for  linear  problems. 

The  proposal  described  12  particular  avenues  of  investigation.  Each 
of  these  will  be  covered  briefly  In  the  next  section. 


SUMMARY  OF  MOST  IMPORTANT 
TECHNICAL  RESULTS 

1.  Fourth-Order  Driver  (FOP)  Methods.  This  type  of  semidirect  method 
was  successfully  developed,  but  not  quite  as  envisioned  in  the  proposal. 

The  intention  was  to  use  the  Bauer-Relss  block-5  code  for  the  linear  solver. 
Instead,  the  marching  methods  were  extended  to  block-5  matrices,  and  an 
even  better  approach  was  developed  solving  the  two  second-order  equations 
for  vortlclty  and  stream  function  In  a coupled  manner.  The  marching  method 
Is  faster  and  allows  much  larger  problems  than  the  Bauer-Relss  method.  The 
resulting  semidirect  nonlinear  method  Is  preferable  to  the  previous  sequential 
equation  approach,  especially  for  strongly  separated  flows,  since  It  requires 
no  semi -empirical  relaxation  factors. 

2.  Methods  for  the  Primitive  Equations.  This  was  accomplished  as  anti- 
cipated, with  no  surprises. 

3.  Methods  for  the  Compressible  Flow  Equations.  This  approach  was  proven 
in  concept,  but  was  not  exercised  In  any  interesting  flow  problems. 

4.  Methods  Based  on  Third-Order  Equations  for  Velocities.  This  approach 


was  superceded  by  the  success  of  #1  above. 


5.  Viscous  1-D  Compressible  Flow  with  Shocks,  This  interesting  problem 
was  unfortunately  not  studied,  due  to  time  limitations  and  due  to  a prior- 
ity of  other  Interesting  and  successful  problem  areas. 

6.  Transonic  Equations.  This  problem  has  been  formulated,  and  we  have 
made  arrangements  with  Dr.  M.  Hafez  to  work  on  the  physical  problem  with 
us.  However,  we  have  postponed  the  final  Implementation  until  we  have  a 
structured  code  for  the  linear  solver  including  stabilization,  so  that  a 
highly  refined  mesh  can  be  used. 

7.  Higher  Order  Solutions.  The  deferred  correction  approach  described 
in  the  proposal  was  tested  and  proven  to  work  quite  well  for  0(a4)  accuracy. 
Attempts  to  achieve  0(a6 ) accuracy  have  failed,  even  In  1-D  linear  problems. 
It  has  also  been  learned  that  gradient  boundary  conditions  should  be  solved 
to  0(a4)  directly,  rather  than  by  the  deferred  correction,  to  enhance  Itera- 
tive convergence.  We  have  not  found  a satisfactory  Implementation  of  a 
high-order  wall  vortlclty  In  the  coupled-equation  approach,  although  it 
does  work  well  In  the  sequential -equation  approach. 

8.  Improved  Boundary  Vorticity  for  Split  NOS  Method.  We  have  success- 


fully Implemented  Israeli's  method  for  wall  vorticity,  and  extended  it  to 
a high  order.  This  method  Improves  accuracy.  Is  easily  applied  in  non- 
orthogonal  coordinate  systems,  and  simplifies  the  evaluation  of  Iterative 
convergence,  so  It  is  now  recoimended  generally  for  the  sequential -equation 
approach.  We  also  experimented  extensively  with  multi-level  iteration 
methods  and  met  with  minor  success  in  Improving  Iterative  convergence,  but 
primarily  we  developed  an  appreciation  for  the  limitations  of  conventional 
one-point  Iteration  schemes  for  wall  vorticity,  which  are  superceded  by 
#1  above  for  difficult  problems. 

9.  Improved  Estimates  of  Advanced  Nonlinear  Terms  by  Quasi linearization 


and  Extrapolation.  We  just  recently  realized  the  limitation  on  any 
quasi linearization  (generalized  Newton)  method  In  regard  to  the  operation 
count.  For  an  MxM  problem,  even  Iterative  Improvement  of  the  Jacobian 
would  require  0(M3)  operations  per  Iteration,  so  that  we  have  settled  on 
the  Split  NOS  or  quasi-Picard  Iteration  for  multi-dimensional  problems. 
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Gartltng  has  recently  performed  a comparison  of  these  iteration  schemes  for 
finite  element  methods  and  has  come  to  the  same  conclusion;  he  also  finds 
the  Split  NOS  method  less  sensitive  than  Newton  iteration  to  initial  guesses, 
especially  at  high  Re.  For  ID  problems,  Newton  iteration  is  clearly  superior. 

Unfortunately,  we  have  not  found  time  to  experiment  with  extrapola- 
tion of  the  nonlinear  coefficients. 

10.  Turbulence  Equations,  Non-Newtonian  Fluids,  and  Temperature-Dependent 

Viscosities.  The  generic  equation  for  all  these  is  a nonlinear  diffu- 
sion term,  typified  by  a turbulence  equation  using  a simple  eddy  viscosity 
formulation.  This  equation  has  been  successfully  implemented.  As  expected, 
the  additional  nonlinearity  caused  slightly  slower  Iterative  convergence 

and  more  sensitivity  to  the  initial  guess,  but  the  method  still  performs  ex- 
cellently. A pleasant  surprise,  easily  analyzed  and  predicted  after  the 
experiment,  is  that  the  addition  of  eddy  viscosity  stabilizes  the  linear 
marching  method.  Also,  a method  for  predicting  optimum  under-relaxation  at 
a wall  was  developed. 

11.  Application  to  Complex  Geometries.  This  has  been  accomplished  with 
considerable  generality.  The  inclusion  of  cross-derivatives  in  the  linear 
marching  methods  has  allowed  general  non-orthogonal  coordinate  transforma- 
tions, and  the  extension  of  Israeli's  boundary  condition  for  wall  vorticity 
to  non-orthogonal  systems  has  allowed  high  accuracy  with  even  some  improve- 
ment in  iterative  convergence.  Codes  have  been  written  and  exercised  which 
allow  arbitrary  tabular  input  of  wall  shape,  with  the  necessary  wall  deri- 
vatives being  obtained  by  0(a6)  differences. 

12.  Three-Dimensional  Problems.  The  success  here  has  been  mixed.  The 
limitations  of  the  direct  marching  method  simply  extended  to  3D  (3D  EVP)  has 
been  defined,  with  a 15x15x30  grid  emerging  as  about  the  upper  limit.  Also, 
the  speed  Is  not  nearly  as  favorable  as  in  2D.  A very  successful  3D  EVP/ 

FFT  linear  solver  has  been  developed,  using  a Discrete  Fourier  Transform  in 
the  3rd  coordinate,  z.  This  unfortunately  restricts  the  equation  solved 

to  terms  like  32/3z2  in  z.  In  the  semidirect  Iteration,  advection  terms 
must  be  lagged.  Under  separate  funding,  it  was  demonstrated  that  the  non- 
linear iterations  using  several  schemes  were  unstable,  apparently  due  to 


j 


I 


3 


J 

! 


the  sensitivity  of  the  FFT  to  small  perturbations.  The  3D  performance  has 
been  the  most  disappointing  aspect  of  the  research. 

As  Is  to  be  expected  In  basic  research,  several  other  developments 
occurred  which  were  not  outlined  in  the  proposal.  A brief  description  of 
these  follow. 

13.  Additional  Developments  in  Marching  Methods  for  Linear  Equations. 

The  basic  tool  used  for  the  nonlinear  semi direct  methods  Is  the  marching 
method  (EVP  method)  for  linear  equations.  Major  developments  were  achieved 
and  described  In  a three-part  review  article  In  Numerical  Heat  Transfer. 

Some  of  these  have  been  described  above;  In  addition,  other  developments 
are  Indicated  here. 

13A.  Accurate  and  complete  operations  counts  were  derived  and  a 30% 
(typical)  savings  was  realized  by  avoiding  some  null  calculations. 

13B.  Stable  methods  for  gradient  and  mixed  boundary  conditions  were 
developed. 

13C.  The  destabilizing  effect  of  upwind  differencing  was  discovered 
and  explained. 

13D.  Cross-derivatives  and  Helmholtz  terms  were  included. 

13E.  A simple  expression  of  previously  graphical  results  for  maximum 
problem  size  was  developed. 

13F.  An  expanding  grid  based  on  the  Fibonacci  sequence  was  developed. 

13G.  Several  rapidly  converging  global -Iterative  methods  for  the 
linear  equations  were  developed,  using  direct  solutions  on  a submesh  and/ 
or  banded  approximations  to  the  full  influence  coefficient  matrix. 

13H.  A general  code  for  the  solution  of  a 9-polnt  or  5-polnt  linear 
operator  was  developed.  It  was  used  by  M.  Luther  of  UNC  on  the  CRAY  compu- 
ter at  NCAR  this  spring,  on  an  ocean  circulation  problem. 

131.  Several  techniques  were  developed  (both  by  the  P.I.  and  by 
other  researchers)  and  compared  for  stabilizing  the  marching  methods.  This 
Is  a major  development  In  the  state  of  the  art,  and  promises  to  make  marching 
methods  Into  robust  tools  for  numerical  solutions  of  physically  important 
problems. 


4 


14.  A peculiar  instability  in  the  Split  NOS  iteration  scheme  was  dis- 
covered. It  arose  for  high  Re  separated  flow  problems  when  the  initial 
guess  on  the  velocity  field  was  too  accurate.  A stability  analysis  showed 
that  the  method  Is  unstable  at  stagnation  points  in  the  initial  velocity 
field;  the  problem  is  simply  avoided  by  requiring  a minimum  value  of  |V°|. 

15.  With  D.  Gartling,  a comparison  of  semi  direct  iteration  schemes  used 
in  finite  difference  and  finite  element  solutions  was  made.  This  Included 
a thorough  operation  count  not  only  for  the  linear  solution,  but  also  for 
the  entire  nonlinear  solution  procedure. 

16.  Vorticity-Velocity  Equation  System.  Like  the  implementation  of  the 
primitive  variable  (pressure-velocity)  equation  system,  this  application 
was  successful  and  held  no  surprises. 

17.  Cylindrical  Coordinate  Problem  with  Exponential  Coordinate  Mapping 
and  Rezone  Mesh  Refinement.  This  code  was  developed  under  partial 

funding  of  this  contract  and  partially  under  a BRL  contract. 

18.  Natural  Convection  Method.  A semi direct  method  has  been  developed 
which  is  especially  suitable  for  natural  convection  problems;  the  stream 
function  and  vorticity  equations  are  solved  in  a coupled  manner,  followed 
by  a (sequentially)  lagged  temperature  equation. 

19.  An  evaluation  was  made  of  the  papers  at  the  First  International 
Conference  on  Numerical  Methods  in  Laminar  and  Turbulent  Flows,  held  at 
University  College  of  Swansea  in  Wales  during  July  1978.  This  evaluation 
focused  primarily  on  work  being  done  overseas,  and  was  written  up  in  a 
trip  report  for  ARO  dated  15  January  1979. 
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